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In this paper, we study a system of two Rossler oscillators coupled through a time-varying link, 
periodically switching between two values. We analyze the system behavior with respect to the 
frequency of the switching. By applying an averaging technique under the hypothesis of a high 
switching frequency, we find that although each value of the coupling does not produce synchro¬ 
nization, switching between the two at a high frequency stabilizes the synchronization manifold. 
However, we also find windows of synchronization below the value predicted by this technique, 
and we develop a master stability function to explain the appearance of these windows. Spectral 
properties of the system are a useful tool for understanding the dynamical properties and the syn¬ 
chronization failure in some intervals of the switching frequency. Numerical and experimental results 
in agreement with the analysis are presented. 


I. INTRODUCTION 

The interaction of two nonlinear (in particular, chaotic) systems via coupling typically leads to a variety of signifi¬ 
cant behaviors, among which the most intriguing is probably synchronization, that is, the coordination of a particular 
dynamical property of their motion [I|. The strength and inevitable delay in the coupling affect the type of coor¬ 
dinated dynamical property, leading to different types of synchronization. The most common forms are complete 
synchronization, i.e. a state in which the two dynamical units evolve following exactly the same trajectory, phase 
synchronization ii when the coordinated property is the phase, lag synchronization Q when both amplitudes and 
phases are locked but with a permanent time lag, and generalized synchronization Q when a given function of the 
output of two systems is synchronized. 

The phenomenon is even more varied when more than two units interact according to a pattern of connectivity 
through which they share information about their current state: complete synchronization , phase synchronization 
0, cluster synchronization Q, partial synchronization @, chimera states [13, HIl: relay synchronization [l3 |. and 
remote synchronization [l3| are only examples of the behaviors observed in a network of coupled oscillators. For 
complete synchronization as investigated in this paper, a series of well-known results provides the conditions for the 
onset of the synchronous motion. If the synchronization manifold is a stable solution for the given network, a suitable 
coupling configuration can always be designed (l^ . 

Although the connectivity among the dynamical units is usually considered time-invariant, interaction among 
dynamical systems may also occur in a discontinuous way (for instance, when it is mediated by links activated 
according to the relative distance of mobile units [l3 - [l0 1 or with weights varying in time according to some adaptation 
law (isl - l^ . In such cases, a key factor in determining the global behavior of the system is the interplay between its 
time scales: one related to dynamics of the units, and the other defining the rate of variation of the links between 
them. 

Another important ingredient is the dynamical rule for the variation of the coupling, which may be either a 
stochastic/periodic activation/deactivation [2l|, H^l or a deterministic law [23|. Both possibilities were explored in 
early experimental works on two coupled Chua’s circuits (^. [^. In particular, in a daptive coupling was used to 
design communication channels able to compensate for parameter changes, while in [25l| it was demonstrated that in 
a synchronization scheme where two Chua’s circuits are pulse coupled, the switching frequency has to be larger than 
a threshold value to achieve synchronization. This behavior is now grounded on recent theoretical results for blinking 
networks [l^, 1131 , proving that under an assumption referred to as the fast switching hypothesis (FSH), when the link 
changes occur enough faster than the oscillator dynamics, the time-varying coupling can be studied by means of the 
time-average of the coupling matrix. 
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However, since there is some arbitrariness in the definition of “/ast enough ,procedures to determine explicit bounds 
for the time-scale of the process driving the coupling mechanism are currently under investigation . On the other 
hand, there is evidence that even below the threshold given by the FSH, many interesting phenomena may occur. For 
instance, a recent study on synchronization of chaotic oscillators coupled via an on-off stochastic network has unveiled 
non-trivial windows of complete synchronization even when, under the FSH, synchronization is not predicted [^ . 

This paper is a case study of two Rossler systems interacting through on-off coupling and showing interesting 
phenomena in a regime not dominated by fast switching. We consider the case where the two chaotic systems are 
coupled through a link whose weight is time-varying, i.e. the weight switches between two fixed values with a given 
switching frequency. We fix the values so that neither of them provides synchronization in the case of static coupling. 
Despite this, switching between the two values gives synchronous behavior. The stability of the synchronization 
manifold is studied for switching frequencies spanning the slow to fast switching regimes, unveiling the relationship 
between the two time-scales. We also provide experimental validation of the results using a hybrid platform where 
two analog electrical circuits are coupled by digitally controlled links. 

The rest of the paper is organized as follows: in Section |IT] the model is described; in Section nil the numerical 
results are shown, and in Section ITVl the analysis with respect to the switching frequency is presented, while Section IVl 
deals with the experimental validation of the results. Section EH gives conclusions. 


II. MODEL 


The Rossler oscillator is described by the following nonlinear dynamical equations 


X = —y — z 

y = x + ay (1) 

z = b + z{x — c) 

where a, b, and c are system parameters. We fix a = 6 = 0.2 and c = 7 throughout the rest of the paper so that a 
chaotic attractor is obtained in the range of considered initial conditions. 

We first consider two identical Rossler systems with a time-invariant diffusive coupling acting between the first 
components of their state space, and we briefly discuss the behavior of this configuration as a prelude to the fast 
switching approach. The equations governing the two coupled systems can be written as 


ii = -yi - zi + k(x 2 - Xi) 

yi = xi + ayi 
zi = b + zi{xi - c) 

X2 = -y2 - Z2 + k{xi -X 2 ) ^ 

y2 = X2+ ay2 
Z 2 = b + Z2{X2 — c) 

where the subscripts indicate the two respective systems, and k is the coupling strength. For brevity, we define 
xi = [ a;i yi zi and X 2 = [ X 2 ?/2 Z 2 ]^, respectively, and rewrite Eqs. m compactly as 


Xi = /(xi) -b kE(x 2 - Xi) 
X 2 = /(X2) -b kE(xi - X2) 


where E = 


1 0 0 
0 0 0 
0 0 0 


Complete synchronization is formally defined as 


||xi — X 2 II —0, as t —>■ 00 (4) 

where || • || stands for the Euclidean norm. To derive the conditions on k for complete synchronization, we define the 
error as e{t) = Xi — X 2 and calculate the error dynamics from Eqs. ([3]): 


= /(xi) - /(X2) - 2 kE(xi - X2) 


( 5 ) 
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FIG. 1: Master Stability Function for the system in Eqs. © with static coupling. 


By linearizing around the common solution xi = X 2 = s, we obtain 


de 

dt 


a/(x) 

9x 


-2ke') 


( 6 ) 


The maximum (largest) Lyapunov exponent Amaxi^) of Eq. ([S|) is a function of k and indicates the region of local 
stability of the error dynamics, and so of complete synchronization. Following the terminology of [l^, we refer to it as 
the Master Stability Function (MSF) of the system. In particular, synchronization requires that Amaxi^) < 0. Figure 
[Udisplays the MSF for the static coupling case of Eqs. ([2]). Note that Ajnax{n) < 0 only in the interval 0.1 < k < 2.35. 
This behavior, referred to as a class-II MSF [s^, is characteristic of a class of systems, including the Rossler oscillator 
when coupled through the variable x. 

We now focus on the main object of the study, which is a system formed by two Rossler units coupled with a 
time-varying link given by 


il = -yi - Zi+ K{t){x2 - Xi) 

2/1 = Xi -I- ayi 
zi = b + zi{xi - c) 

±2 = -2/2 -22-1- K{t){xi -X 2 ) ^ 

y 2 =X 2 + 02/2 

Z 2 = b + Z2{X2 - c) 

We assume the coupling is given by K{t) = ki -f (sgn(cos(a;t)) -|-1) where sgn(a;) = 1 for a; > 0 and sgn(a:) = —1 
otherwise, so that the effective coupling switches between two constant values fci and k 2 at a frequency w. We refer 
to Eqs. o as the switching system and analyze its behavior with respect to the switching frequency uj, which is an 
important bifurcation parameter. 

In particular, we select ki and ^2 such that neither of the two falls within the synchronization range for the static 
coupled case (MSF of Fig. [1]), that is, Amax{ki) > 0 and Amax{k 2 ) > 0. Under these conditions, the problem of 
synchronization of the switching system is not trivial since the system switches between two configurations that are 
not synchronizable. A recently developed approach for blinking systems provides a useful tool for understanding 
the behavior of the system when the switching occurs at a high enough frequency. For an average system defined 
by substituting for the coupling parameter its average value (in our case k{t) = k = ), the trajectory of the 

switching system approaches that of the average system under the hypothesis of fast switching. We will show that 
the fast switching approach can be used to predict the behavior of the switching system, although synchronization 
can also occur at lower frequencies within windows of uj. 


III. NUMERICAL RESULTS 

In the numerical simulations, we set fci = 0 and k 2 = 2.4. Both coupling values fall outside the range of synchro¬ 
nization identified by the MSF_of Fig. [H i.e., Amax{ki) > 0 and Amax{k 2 ) > 0, while the average value k = 1.2 lies in 
the stable region where Amax{k) < 0. 
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FIG. 2: Average synchronization error E{ijj) with respect to the switching frequency uj. 




FIG. 3: Complete synchronization of two Rossler oscillators coupled by a switching signal with u) = 1.5. Projection of the 
attractor onto the planes (a) xi-y\ , (b) X2-y2, (c) Xi-X2- 


We investigate the effect of the switching frequency by fixing all the other parameters and varying uj from 0 to 1.5 
(the limit case a; = 0 corresponds to uncoupled dynamics). For each value of w, we integrate Eqs. ([7]) with a 4th order 
adaptive Runge-Kutta algorithm for a time T = lO^s and sample the result at dt = 0.01, thus obtaining M = 10® 
samples from which the average synchronization error E(lo) is calculated from 


\/{xi{h) - X2{h)f + {yi{h) - y 2 {h)Y + (zi(fe) - Z2{h)f 

yJxlihY + yiQiY + Zi{hY + X2{hY + y2{hy + Z2{hY 

The synchronization error E{ui) is normalized so that E = 1 means that the two trajectories are completely 
uncorrelated, E > 1 indicates anti-correlation, while E ^ 0 corresponds to the highest correlation. 

The synchronization error E{u!) in Fig. [5] is a nonmonotonic function of uj. In particular, for w > 1.3 the system 
is synchronized. Therefore, w ~ 1.3 represents a boundary between the region in which the FSH holds, and that in 
which the frequency of the switching is not “fast enough.” For uj > 1.3, the time-varying connectivity is sufficiently 
fast relative to the Rossler dynamics that the system responds as it would to a constant coupling equal to the average 
of the two coupling strengths fci and ^ 2 - In fact, the two oscillators in this regime are always synchronized and chaotic 
as shown in Fig. [3] for w = 1.5. 

The behavior at low frequencies {uj < 0.2) is expected from the choice of fci and k2. In this case, the system slowly 
alternates between two configurations, neither of which is synchronizable. The global behavior is thus unsynchronized. 

The most interesting frequencies are in the range uj € [0.2,1.3] where nontrivial alternating windows of synchro¬ 
nization and nonsynchronization occur. In particular, even when the FSH does not hold and the oscillating link has 
time scales comparable to those of the Rossler system, there are regions of complete synchronization. 

Starting from a; = 1.3 at the boundary where the FSH holds and decreasing the switching frequency, we examine 
the dynamical behavior of the Rossler oscillators at different values of uj. Just below ui = 1.3, a large window 
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FIG. 4: Unsynchronized behavior of two Rossler oscillators coupled by a switching signal with ui = 1. Projection of the attractor 
onto the planes (a) xi-yi, (b) X2-y2, (c) Xi-X2. 



FIG. 5: Temporal evolution of the average synchronization error E for lj = 1.13 showing intermittency. 


(0.9 < u! < 1.3) of unsynchronized behavior is found. The system attractor is significantly different from the chaotic 
attractor of the synchronized Rossler oscillators as shown in Fig. |4] for uj = 1.0. 

We also observe two other significant dynamical behaviors. At w 0.97 and w w 1.13, which correspond to the 
transition between synchronous and unsynchronous regions, intermittent synchronization occurs as shown in Fig. [S] 
The shape of the attractor in this case alternates between the original shape of the Rossler chaotic attractor (obtained 
when the synchronization error is close to zero) and the one shown in Fig. |4l obtained when the error is larger. 
Furthermore, at a; ~ 0.9656 there is a narrow window in which two stable limit cycles coexist as shown in Fig[51 

For a further decrease in the switching frequency, a new window of synchronization occurs at 0.6 < uj < 0.9. Below 
w = 0.6 a series of unsynchronized/synchronized windows are observed whose widths decrease for decreasing values 
of UJ. Within all these windows, including the main one around w = 1, a narrow window of synchronization occurs, 
only three of which are evident in Fig. [51 

The observed behavior is confirmed by an analysis of the four largest Lyapunov exponents ([7]) shown in Fig. [71 
Regions of synchronization are characterized by one positive Lyapunov exponent, while those of unsynchronized 
behavior have two positive Lyapunov exponents corresponding to hyperchaos in the seven-dimensional state space. 
The sign of the second largest non-zero Lyapunov exponent thus discriminates between the windows of synchronous 
and unsynchronous motion found by analysis of the synchronization error. 


IV. ANALYSIS 

From the numerical results of Section uni it is evident that the switching frequency is a bifurcation parameter, 
especially when the FSH does not hold, that is, when the rate of switching is comparable to the dynamics of the 
Rossler system. The switching frequency affects not only the synchronization but also the dynamics of the attractor. 
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FIG. 6: Limit cycles in the two Rossler oscillators coupled with a switching signal with uj = 0.9656. (a) Projection of the 
attractor onto the plane xi-yi, (b) projection of the attractor onto the plane X 2 -y 2 , and (c) projection of the attractor onto the 
plane Xi-X2- 



FIG. 7: Lyapunov exponent spectrum for the system in Eqs. 0. 


This section provides further analysis of the observed phenomena. In particular, we first investigate the properties of 
the power spectra for the Rossler state variables and correlate them with the windows observed in Figs. [2] and [7] and 
then develop a MSF for the case of two oscillators with a time-varying coupling. 

We begin by observing that the power spectral density of an uncoupled Rossler system 0 as shown in Fig. [5] is 
characterized by a strong dominant component corresponding to the large nearly periodic oscillations in x-y near the 
z = 0 plane having a frequency ujs — 1.067. The large window of unsynchronized behavior in Fig. [2] just below the 
boundary with the fast switching region is around this frequency. A closer inspection of the switching system spectrum 
shows that a resonance occurs in this window as shown in Fig. [S) The Figure shows the power spectral density (color 
coded) for 0.9 < lUs < 1.2 as a function of the switching frequency lu € [0.5,1.5]. In the two windows where the 
systems are synchronized (one is the fast switching region for w > 1.3, and the other is the window 0.6 < w < 0.9), the 
dominant frequency is ujh, the same as for the uncoupled Rossler system. That is, when the two systems synchronize, 
they evolve following the chaotic trajectory s of the uncoupled system. On the contrary, when the coupling strength 
switches at a frequency oj comparable to ojji, in particular in the window 0.95 < w < 1.15, the dominant frequency 
is locked to the switching frequency, thus resulting in a modification of the dominant frequency with respect to the 
uncoupled case. This explains the different shape of the attractor in this range (Fig. 0]). A similar locking occurs at 
the first subharmonic of the dominant frequency around oj = 0.5335 superimposed on the window of unsynchronized 
behavior. Presumably, the same locking occurs in narrower ranges around all the other subharmonics. 

The frequencies of the three unstable periodic orbits with the shortest periods for the single Rossler system are 
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FIG. 8: Power spectral density for the x variable of the Rossler system in Eq O- 
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FIG. 9: Power spectral density for the x variable of the Rossler system in Eq (l7|), for different values of tu. In the color code 
red indicates higher power. 


1.0626,0.5294, and 0.3533, corresponding closely to the dominant frequency and its first two subharmonics. Thus it is 
difficult to distinguish whether the resonances are a result of the peak in the power spectral density or an interaction 
with the unstable orbits on the attractor. However, it may be significant that the number of unstable periodic orbits 
proliferate enormously for frequencies below about w = 0.3, and this may account for the absence of low-frequency 
structure in E{oj) shown in Fig. [2j 

We now derive a MSF for the time-varying coupling. To do this, we calculate the error dynamics for system dZD 
and linearize around the common solution Xi = X 2 = s to obtain 


d{e) 

dt 


( a/(x) 
V ax 




(9) 


with K{t) = ki + (sgn(cos(a;t)) -b 1). From Eqs. ([5]), the maximum Lyapunov exponent Amaxi^j) as a function 
of the parameter oj is calculated as shown in Fig. [TUI Values of w such that Ajnaxi^) < 0 correspond to switching 
signals able to synchronize the Rossler oscillators, while A^axioj) > 0 indicates that the error does not decay to zero 
















FIG. 10: Maximum Lyapunov exponent Kmaxiui) for the system in Eqs. 0 with time-varying coupling. 


for that value of the switching frequency. The curve Amaxi^j) vs. w perfectly explains the windows of synchronization 
and nonsynchronization for system 0. 


V. EXPERIMENTAL INVESTIGATION 

In this section, the system of two Rossler oscillators with switching coupling is experimentally investigated using 
an electronic circuit governed by Eqs. 0. Each Rossler oscillator uses the electrical scheme reported in [3l|. The 
temporal dynamics are rescaled by a factor K = 2128, so that the circuit waveforms correspond to those of the Rossler 
system o with a rescaled time axis. 

Unlike static coupling, which requires only a single resistor between the capacitors associated with the respective 
coupled state variables, coupling in this experiment used an analog switch in series with a coupling resistor driven by 
a pulse width modulated (PWM) source. This strategy is implemented in anticipation of extending the studies to a 
network of switched connections. We have also explored solutions using a digital resistor, but the resolution provided 
by off-the-shelf digital resistors in terms of values and switching time is inadequate for our purposes. 

In more detail, the coupling scheme was implemented by using two components: the analog switch ADG452 and 
the ST microcontroller unit (MCU) STM32F303VCT6 for the generation and control of the PWM signal. The 
ADG452 embeds four independently selectable bi-directional switches, has a low on-resistance (on the order of 517), fast 
switching times {toN = 70ns, toFF = 00ns), and is TTL-/GMOS-compatible. The microcontroller STM32F303VGT6 
is an ARM®-based Gortex@-M4, 32bit microcontroller with an embedded floating point unit. It has a core clock 
speed up to 72MHz, a 256KB flash memory, 48KB SRAM, and a wide range of peripherals such as analog-to-digital 
and digital-to-analog converters, timers, direct memory access, etc. It operates with a voltage supply in the range 
from 2.0U to 3.6U. 

The value of the coupling resistor is controlled by the duty cycle (DC) of the PWM signal driving the analog switch. 
Turning on and off the switch has the effect of multiplying the fixed coupling resistor by a factor inversely proportional 
to the DC, according to the equation 





(10) 
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FIG. 11: Average normalized error E{uj) calculated on data acquired from the experimental system. Frequency axis is rescaled 
by the factor K to allow comparison with numerical results. 

In our application, the ADG452 has been controlled with a 40kHz PWM, which is a suitable value since the 
characteristic frequency range of the Rossler implementation is below 5kHz. Although this solution is able to switch 
between two nonzero values of the coupling resistance, it allows implementing a very low value of the coupling ki 
{ki = 0.02), which for the purpose of our analysis is equivalent to two disconnected circuits. 

The waveforms corresponding to the six state variables generated in the experiment were acquired using an NI- 
USB6255 data acquisition board at a sampling frequency of fs = IbkHz and post-processed to compute the average 
normalized error E{fl) as in Eq. ([5]), where Q = Koj is the switching frequency in the rescaled circuit. 

The trend of E{n) is shown in Fig.[TTJ To facilitate the comparison with the numerical results, the frequency has 
been rescaled in terms of the variable uj. The presence of a main peak of desynchronization around lo = 1.06 is also 
confirmed in the experimental case. Due to component tolerances, in spite of predicted complete synchronization, i.e., 
E = 0, only practical synchronization (sH, {E < 0.2V) is observed. The windows of practical synchronization are 
in good agreement with those of complete synchronization found in the numerical simulations. 


VI. CONCLUSIONS 


In this paper, two Rossler oscillators coupled through a time-varying link have been investigated. In particular, the 
coupling strength switches periodically between two values. The observed behavior is significantly more rich than the 
case of other chaotic circuits such as Chua’s circuit, for which there is only a single transition in frequency between 
synchronized and unsynchronized behavior 


25|. 


Although each of the two values of the coupling, when applied statically to the system, does not lead to synchro¬ 
nization, switching between them can stabilize the synchronization manifold. This behavior is observed not only for 
high switching frequencies as predicted by the fast switching theory, but also in several windows at lower switching 
frequencies. 

The alternate windows of synchronized and unsynchronized behavior can be explained by a MSF illustrating the 
behavior of the maximum Lyapunov exponent transverse to the synchronization manifold as a function of the switching 
frequency. Windows of synchronization correspond to negative values of the maximum transverse Lyapunov exponent, 
while unsynchronized behavior is obtained when this measure is positive. 

We have also demonstrated a strong effect of the switching coupling on the spectral properties of the two oscillators, 
since a switching frequency very close to the dominant component of the spectrum of the isolated Rossler oscillators 
causes the dominant component to lock to the switching frequency, resulting in a significant change of the chaotic 
dynamics and a failure of synchronization. 

Hence inspection of the spectral properties of the system and analysis of the MSF derived for the time-varying 
coupling allow one to establish the regions of synchronizability beyond that predicted by the fast switching analysis. 
These regions agree well with the results of numerical analysis as well as experimental results. 
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